Theory of alike selectivity in biological channels 
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We introduce a statistical mechanical model of the selectivity filter that accounts 
for the interaction between ions within the channel and derive Eisenman equation 
of the filter selectivity directly from the condition of barrier-less conduction. 


Selectivity of the canonical K* filter, formed by a narrow 15 A-long tunnel holding 
four Kt-binding sites, is usually analyzed [1-6] in terms of the difference in the excess 
chemical potential in the bulk /; in and the channel jf at one binding site 


AAjix,Na = (Hina — Fina) — (xk — HK): (1) 
The question of how to introduce barrier-less conduction within this formalism 


elicited a many-voiced debate over decades. 


Let n? be a number ions of the 7-th type in the solution and n,; their number within 
the filter. We introduce the free energy of the solution plus the filter in the form 


G = So(ni — n)kT naj + SO(n} — nu) fii + DO nifif + kT In (11 ni) +e({n;},n-). 


In treating the filter filled with ions as excited states of the system, we derive the 
following grand canonical ensemble for the filter (cf [?]) 
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The current with a constant voltage drop Ad = ¢” — ¢* through the channel in 
linear response is proportional to the rate of charge fluctuations passing through the 


channel [7] J steady = ee Ad. 


For discrete conducting ions, charge fluctuations are proportional to fluctuations 
in the number of ions in the channel. The variance 02, = —kT es of the 
t/ TV 


number of ions in the filter can be found using the grand potential Q = —kTInZ. 


To apply obtained results to the KcsA filter, let us consider two types of conducting 
ions Kt and Na* that can bind to any binding site in the filter, but have to have at 


least one water molecule between them. The free energy of the filter filled with nx 
and nwa ions (ignoring a prefactor) is G = —nK Ajix — nnaAfina + €(Nk, Na, NF). 


The energy of interaction between the ions is e(nx,NNa, Nf) = r (nk +NnatMNf), 


where n,- is the number of fixed charges on the channel wall, i.e. the main mutation 
parameter. The channel capacitance is C & Ament Re 

The lowest conducting free energy level corresponds to nx = 2. The conditions for 
fast conduction are: (i) G(ng,n+) = G(ng + 1,n}) and (ii) G(ng,n}) © OkT. To 
tune the selectivity and conductivity of the filter one has to tune: (i) geometry, (ii) 
jig, and (iii) the ns. Tuning the filter for fast conduction of K* ions gives 


C 
ny = —(nxe + 1/2)+ GAue. (3) 


Substituting this value into the free energy gives the following barrier for Nat ion to 
enter the channel AGy, = Aux — Apna, i.e. Eisenman’s selectivity relation follows 
directly from the condition for fast conduction of K* ions. 


Analysis of the mean and variance for Kt and Na* ions in the filter shows that 
tuning it for fast conduction of Kt ions results in effective block of Na* conduction 
if the width of the dispersion peaks is smaller than the separation of the peaks 
determined by the geometry of the channel and AG. 


In conclusion, we have introduced a statistical mechanical model for the distribution 
of alike-charged ions in the K* selectivity filter and have shown that the Eisenmann 
relation for filter selectivity follows directly from the condition for fast barrier-less 
conduction of Kt ions in the filter. The results can be extended beyond linear 
response using a solution of the Nernst-Planck equation. Models of filters with 
distinguishable binding sites will be considered in future work. 
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